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摘 要 : 运筹 学 研究 领域 中 的 应 急 服务 设施 选 址 问题 有 许多 求解 模型 ， 该 文选 取 了 P- 中 心 模型 进行 研究 ， 首 先 研究 
了 该 问题 的 数学 性 质 ， 并 给 出 证 明 ， 利 用 这 些 数 学 性 质 能 对 问题 进行 降 阶 从 而 缩小 问题 的 规模 ， 然 后 在 此 基础 上 设 
计 一 个 基于 上 界 和 下 界 的 回溯 算法 来 求解 该 问题 ， 最 后 通过 一 个 示例 分 析 进 一 步 阅 述 该 算法 的 原理 ， 并 证 明了 该 算 
法 能 在 较 短 时 间 内 求 得 问题 的 最 优 解 。 
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Backtracking algorithm with reduction for p-center location problem 


Shang Chunjian, Ning Aibing, Peng Dajiang, Zhang Huizhen 
(Business School, University of Shanghai for Science & Technology, Shanghai 200093, China) 


Abstract: There are various solution models for the emergency service facility location problem and this paper studies the P- 
center model. This paper proposed a backtracking algorithm with upper and lower bound, and the algorithm added a reduction 
algorithm by studying the mathematical properties. The algorithm can decrease the scale and the degree of complexity of the 
problem, so as to accelerate the execution speed. At the end of the paper, this paper illustrated an instance to elaborate this 
algorithm further. 
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0 ”引言 学 者 研究 精确 算法 求解 P- 中 心 问 题 ， 例 如 Sandoval 等 人 090 

而 提出 基于 整数 规划 模型 的 精确 算法 ，Contardo 等 00 引 入 了 一 
日 常生 活 中 的 各 类 紧急 突 发 事件 给 人 们 的 安全 和 生活 造 ”种 可 伸缩 的 基于 松弛 的 迭代 算法 设计 一 种 精确 算法 求解 P- 
成 了 巨大 的 影响 ， 当 灾害 发 生 时 ， 应 急 服务 设施 是 提供 救援 。 ”中 心 问题 。 
支持 、 保 障 人 民 群 众人 身 安全 最 基本 最 重要 的 设施 ， 具 有 十 该 问题 的 启发 式 算法 和 近似 算法 的 优点 是 能 够 快速 的 得 
分 重要 的 作用 。 到 可 行 解 ， 不 足 之 处 在 于 容易 陷入 局 部 最 优 02. 31， 该 问题 的 
多 数 应 急 设施 选 址 问题 是 NP 难题 ， 有 多 种 不 同 的 模型 ” 精确 算法 优点 是 能 得 到 最 优 解 ， 不 足 之 处 在 于 不 适合 求解 大 
来 解决 这 类 问题 。 针 对 应 急 设 施 选 址 问题 的 公平 性 要 求 ， 本 ”规模 问题 ， 具 有 较 高 的 时 间 复 杂 度 中。 针对 上 有 算法 的 不 足 
文选 取 P- 中 心 模型 中 进行 研究 ， 它 是 离散 位 置 模 型 的 一 种 特 。 之 处 ， 本 文 设 计 了 一 个 精确 算法 ， 首 先 研究 P- 中 心 模型 的 数 
殊 类 型 ， 是 指 从 nn 个 备 选 的 设施 中 最 多 选取 p 个 设施 ， 对 于 ”学 性 质 ， 并 给 出 相应 的 数学 证 明 ， 利 用 这 些 数 学 性 质 能 批量 
每 一 个 应 急 点 ei 均 由 选中 的 设施 中 距离 其 最 近 的 设施 服务 ， 某 些 确定 开设 或 不 开设 的 设施 ， 进 而 缩小 问题 的 规模 ， 最 后 
要 求 所 有 应 急 点 与 对 应 为 其 提供 服务 的 设施 之 间 的 最 大 服务 。 设计 一 个 基于 上 界 和 下 界 的 回溯 算法 ， 求 得 问题 最 优 解 的 同 
距离 最 小 ， 即 尽量 减少 所 有 应 急 点 与 它们 的 服务 设施 之 间 的 。 ”时 能 利用 数学 性 质 加 快 问题 求解 速度 。 
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最 大 距离 。 在 这 个 模型 中 ， 设 施 没有 容量 约束 。 目 前 有 许多 1 ”问题 描述 

学 者 致力 于 此 选 址 模型 的 研究 ， 有 的 学 者 结合 图 论 知 识 对 问 

题 模型 进行 求解 ， 例 如 HockbaumD) 提 出 一 个 时 间 复 杂 度 为 本 文 用 二 分 图 G = (E, 内 来 表示 该 问题 ,已 = {ei|i=1， 
O(ElogE) 的 图 的 顶点 多 中 心间 题 的 2-opt 算法 , Dyer8 提 出 一 。 2,…,m} 表 示 所 有 应 急 点 的 集合 ， 用 ei 代表 第 i 个 应 急 点 ，m 
个 时 间 复 杂 度 为 O(np) 的 2-opt 启发 式 算法 ; 还 有 的 学 者 对 于 ”代表 应 急 点 的 总 个 数 ; = {1 = 1, 2, …, n} 表 示 设 施 的 集 
P- 中 心 问题 采用 了 模糊 集 的 理论 进行 求解 ， 例 如 Wen 中 构建 ” 合 ， 用 代表 第 j 个 设施 ，n 代表 设施 的 个 数 。 

设施 选 址 模型 时 考虑 到 模糊 需求 ， 实 现 应 急 设 施 选 址 的 完全 本 文 引 入 一 个 额外 的 决策 变量 2 来 表示 所 有 应 急 点 和 相 


覆盖 ; 有 的 学 者 采用 启发 式 搜索 的 方法 求解 , 例如 Domschke ”应 为 其 提供 服务 的 设施 之 间 的 最 大 距离 。 己 知 从 设施 中 最 多 
等 人 四] 采用 二 分 法 的 思想 求解 出 中 心 值 等 于 p 的 + 履 盖 问题 ， ”选取 p 个 设施 ,同时 给 定 每 个 设施 与 应 急 点 的 相应 位 置 关系 ， 
朱 燕 等 {9 结合 层次 分 析 法 ， 采 用 贪 禁 取 走 启发 式 算法 对 p- 中 且 p<n; dj 表示 第 i 个 应 急 点 ei 到 第 j 个 设施 有 的 距离 ， 设 
心 模型 进行 求解 ， 应 用 于 航空 救援 服务 ;Davidovi 等 中] 采用 ”决策 变量 ;: 
蜂 群 算法 求解 P- 中 心 问题 ; 许 户 宸 等 由 利用 改进 KK-Means 算 
法 和 重心 法 相 结合 的 方式 对 P- 中 值 问题 进行 求解 ， 黄 书 强 等 | 
外 提出 一 种 加 权 距 离 连 续 K 中 心 选 址 问题 求解 方法 ， 还 有 些 


闻 


1 应急 点 @ 由 设施 提供 服务 
0, 应 急 点 @ 不 由 设施 提供 服务 
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1 设施 f; 被 选取 
ee 
7 “10, 设施 /不 被 选取 


则 P- 中 心 选 址 问题 具体 模型 如 下 : 


min 0 

yj = 上 i=1,2,.....m (D 

jel jlfieF} 
yij — Xi < 0, RE ep /i Ie n(2) 
xX; < p, G3) 

S.t. 4 jeljljer} 
dy;; -QO<0, i=1,2,.....m (4) 

jel jljier} 
xX) e{0,1}, je 2 (5) 
ye{0d), lm j=12,....n(6) 


目标 函数 2 表明 要 最 小 化 每 一 个 需求 点 与 它 的 服务 设施 
之 间 的 最 大 距离 。 
约束 (1) 表 示 有 且 仅 有 一 个 yij 是 非 零 的 ， 即 有 且 仅 有 
个 设施 可 以 为 一 个 应 急 点 服务 ,保证 了 任意 应 急 点 ei 都 能 找 
到 一 个 设施 为 其 服务 ， 满 足 完 全 覆盖 的 要 求 。 
约束 (2) 表 示 当 部 =0 必 有 J= 0， 即 只 有 被 选中 开设 的 
设施 才能 为 应 急 点 服务 。 
约束 (3) 表 示 最 多 开设 p 个 设施 。 
约束 (4) 表 示 0 必须 大 于 第 i 个 应 急 点 与 它 的 服务 设施 之 
间 的 距离 ， 由 于 对 每 个 应 急 点 均 有 这 个 约束 ， 即 2 必须 大 于 
任 一 应 急 点 与 其 服务 设施 的 距离 。 
约束 (3) 和 (6) 表 示 变 量 的 取 值 为 0 或 1。 
该 问题 是 NP-Hard 问题 ， 除 非 P=NP， 和 否则 该 问题 不 存 
在 多 项 式 时 间 的 精确 算法 pa。 
2 ”数学 性 质 及 子 算法 设计 
2.1 数学 符号 
为 了 便于 描述 ， 本 文采 用 以 下 数学 符号 来 描述 : 
ei，f: 表示 第 i 个 应 急 点 ; 有 表示 第 j 个 设施 ; 
m，n: 表示 应 急 点 的 个 数 ，n: 表示 设施 的 个 数 ; 
= {ei|i= 1,2,…, m}: 表示 所 有 应 急 点 的 集合 ; 
下 = {f17=1,2,…,n}: 表示 设施 的 集合 ; 
p: 表示 最 多 选取 p 个 设施 (p < nn); 
dij: 表示 第 i 个 应 急 点 ei 到 第 j 个 设施 的 距离 ; 
mini_di(i): 表示 应 急 点 ei 到 所 有 设施 的 最 小 距离 ; 
mini_n1(i): 表示 所 有 设施 中 ,设施 mini_m1(i) 到 应 急 点 ei 
的 距离 最 小 ; 
mini_d2(i): 表示 应 急 点 e 到 所 有 设施 的 次 小 距离 ; 
mini_n2(i): 表示 所 有 设施 中 ,设施 mini_n2(i) 到 应 急 点 ei 
的 距离 次 小 ; 
maxi_qd1i( 让 ， 表 示 应 急 点 ei 到 所 有 设施 的 最 大 距离 ; 
minj_d10): 表示 设施 到 所 有 的 应 急 点 的 最 小 距离 ; 
maxj_d1( 四 ， 表示 设施 万 到 所 有 的 应 急 点 的 最 大 距离 ; 
min_di(i, 他 : 表示 应 急 点 ei 到 集合 中 所 有 设施 中 的 最 
小 距离 ; 
min_mi(i, 他 ): 表示 在 集合 对 的 所 有 设施 中 ,设施 min_ni(i 
加 到 应 急 点 ei 的 距离 最 小 ; 
min_d2(i, 他 : 表示 应 急 点 ei 到 集合 全 中 所 有 设施 中 的 次 
小 距离 ; 
min_n2(i, 了 ): 表示 在 集合 对 的 所 有 设施 中 ,设施 min_n2(i, 
习 到 应 急 点 ei 的 距离 次 小 ; 
1: 一 定 开设 的 设施 集合 , 若 FF 中 任 一 设施 不 开设 则 不 
E 取 得 最 优 解 ， 为 全 局 变量 ; 
Fo: 一 定 不 开设 的 设施 集合 , 若 Fo 中 任 一 设施 开设 则 不 
能 取得 最 优 解 ， 为 全 局 变量 ; 
三 =F/ Fi/ Fo: 暂时 未 确定 是 否 开设 的 设施 集合 ， 为 全 


Pal 


施放 入 集合 FFo 中 ， 为 全 局 变量 ; 
FFs= Fs/FF1/ FFo: 回 济 算法 分 情况 时 暂时 未 处 理 的 设 
施放 入 集合 FFs 中 ， 为 全 局 变量 ; 


为 全 局 变量 ; 
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FF 本 中 的 设施 在 回溯 算法 分 情况 时 假设 开设 的 设施 
放 入 集合 FF 中， 为 全 局 变量 ; 


FFo: Fs 中 的 设施 在 回溯 算法 分 情况 时 假设 不 开设 的 设 


best_q: 回调 算法 在 当前 状态 下 已 知 最 好 的 目标 函数 


开设 设施 集合 ， [Fpesd <p, 为 全 局 变量 ; 


值 入 


Fpest: 算法 在 当前 状态 下 已 知 最 好 的 目标 函数 值 对 应 的 


u_q: 算法 在 某 一 状态 下 的 目标 函数 值 上 界 ， 即 1 和 


中 设施 开设 时 的 上 界 ， 为 局 部 变量 ; 


好 it 


o 


2.2 


cur_n: 当前 累计 开设 设施 数 ， 为 局 部 变量 ; 
cur_i: 回溯 算法 的 当前 搜索 层 数 ， 为 局 部 变量 ; 


opt flag: 若 opt fiag=1, 说 明 当 前 算法 已 取得 最 优 解 ; 
opt_flag=0， 不 能 确定 当前 是 否 已 取得 最 优 解 ， 为 全 局 变 


数学 性 质 


J 


性 质 1 车 存在 一 设施 f, 假如 设施 万 开设 且 此 时 的 上 界 


u_g 小 于 当前 的 下 界 best_q, 则 设施 万 一 定 不 开设 , 其 中 上 下 


界 求解 算法 见 后 面 2 节 的 内 容 。 
证 明 上 界 u_g 小 于 当前 的 下 界 best g， 说 明 此 时 开设 
设施 不 可 能 获得 最 优 解 ， 所 以 应 该 关闭 设施 f。 
性 质 2 车 存在 某 一 设施 f, 假如 设施 不 开设 且 此 时 上 
界 u_g 小 于 当前 的 下 界 best g， 则 设施 一定 开设 。 


设施 有 不 可 能 获得 最 优 解 ， 所 以 应 该 开设 设施 ff。 
性 质 3 当 p= 1 时 ， 则 此 时 开设 maxj_q10) 最 小 的 设施 
即 可 。 


证 明 上 界 wu_g 小 于 当前 的 下 界 best 9g， 说 明 此 时 关闭 


证 明 当 p=1 时 , 即 此 时 只 要 选取 1 个 设施 , 选 maxj_d10) 
最 小 的 设施 能 使 目标 函数 最 小 。 


性 质 4 存在 设施 ， 有， 若 对 于 任意 应 急 点 ei € 已 ， 


= 


dij< din, 则 设施 所 一 定 不 开设 , 即将 设施 访 放 入 集合 Fo 中。 
证 明 因为 此 时 连接 到 设施 轧 的 应 急 点 都 可 以 连接 到 设 


施 太 上 且 目 标 函 数 相 同 或 更 好 。 


的 设施 都 不 影响 目标 函数 值 。 


re 
于 


性 质 5 存在 设施 万 若 其 minj_q10) 大 于 等 于 其 他 任 一 


设施 万 的 maxj_di(h)， 则 设施 一 定 不 开设 ， 即 将 设施 万 放 
入 集合 Fo 中 


证 明 因为 此 时 连接 到 设施 访 的 应 急 点 都 可 以 连接 到 设 


施 太 上 且 目 标 函 数 相 同 或 更 好 。 


性 质 6 存在 两 个 应 急 点 er er, 若 对 于 任意 设施 feF， 


有 dij> dy， 则 可 先 删 除 应 急 点 ex, 算法 最 后 再 把 应 急 点 ei; 连 
接 到 任 一 开设 的 设施 都 不 影响 目标 函数 值 。 


证 明 ”由 于 对 于 任意 设施 feFf， 有 dij> dr;， 因 此 把 应 
急 点 ea 连接 到 任 一 开设 的 设施 都 不 影响 目标 函数 值 。 
性 质 7 存在 两 个 应 急 点 ep ert; 车 mini_di(i)> maxi_di( 有 D)， 
则 可 先 删除 节点 et， 算 法 最 后 再 把 应 急 点 ex 连接 到 任 一 开设 


证 明 由 于 必须 服务 应 急 点 e， 因 此 目标 函数 一 定 大 于 


于 mazzz_di(D， 而 应 急 点 ex 到 所 有 设施 的 最 大 距离 小 于 等 


于 mini_d1i(i)， 即 应 急 点 er 到 所 有 设施 的 最 大 距离 一 定 小 于 


等 于 目标 函数 ， 因 此 把 应 急 点 et 连接 到 任 一 开设 的 设施 都 不 


向 目标 函数 值 。 


性 质 8 存在 设施 4， 车 其 maxj_qd10) 小 于 等 于 其 他 任 一 


个 设施 广 的 minj_qi(h)， 则 设施 必定 开设 ， 将 其 放 入 全 


下 中 ， 且 最 优 函 数值 即 为 maxj_d10)。 
证 明 此 时 所 有 应 急 点 都 连接 到 设施 f 即 可 取得 最 优 
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生 质 9 ” 若 当 前 未 被 服务 的 应 急 点 个 数 小 于 等 于 p-| 下 |- 
IFF1|， 此 时 对 于 每 一 个 未 服务 应 急 点 we， 把 应 急 点 ei 连接 到 
设施 min_ni(i, 了 -Fo- FF0) 即 可 ， 这 样 处 理 当 前 未 服务 应 急 点 
是 最 优 的 选择 。 
证 明 ”这样 操作 开设 的 设施 总 个 数 小 于 等 于 p 且 当 前 未 
服务 应 急 点 都 是 连接 到 距离 其 最 近 的 设施 上 ， 这 样 处 理 当 前 
未 服务 应 急 点 是 最 优 的 选择 。 
性 质 10 ”将 每 一 个 应 急 点 i 都 与 距离 其 最 近 的 设施 相连 
接 , 若 此 时 开设 的 设施 总 数量 小 于 等 于 p, 则 此 时 求 得 最 优 解 。 
证 明 显然 成 立 。 
性 质 11 将 应 急 点 mini_di(i)) 值 从 大 到 小 排序 ， 然 后 依 
次 将 应 急 点 ei; 连接 到 设施 mini_d1(i)， 直 到 开设 的 设施 数 为 p 
或 者 所 有 应 急 点 都 被 服务 为 止 ， 之 后 若 还 有 应 急 点 没有 被 服 
务 ， 则 把 这 些 应 急 点 连接 到 已 经 开设 的 p 个 设施 中 距离 该 应 
急 点 最 近 的 设施 上 ， 若 这 些 应 急 点 连接 的 距离 都 小 于 或 等 于 
则 此 时 取得 最 优 解 ，opt flag = 1。 
标 函 数值 为 下 界 max(mini_q1(?)， 因 此 成 立 。 


[ 茎 


芭 


max(mini_d1(?)) 

证 明 此 时 
2.3 下 界 子 算法 

本 文 设计 的 下 界 子 算 法 具体 步骤 如 下 ; 

a) 初始 化 ， 对 任意 ei EE，f EF， 让 每 一 个 应 急 点 均 
所 有 设施 中 距离 其 最 近 的 设施 服务 ; 

b) 判断 当前 设施 开设 个 数 。 若 当前 设施 开设 个 数 小 于 等 
于 p， 此 时 可 求 得 最 优 解 ， 最 优 函 数值 CO=max(mzzi_di(D))(e' 
E EB)， 输 出 目标 函数 值 ， 算 法 结束 ; 若 此 时 开设 设施 个 数 大 
于 p， 此 时 令 max(mini_qdi(i)) = bo， 执 行 下 一 步 ; 

c) 令 了 为 前 面 开 设 的 所 有 设施 的 集合 , 了 中 设施 的 个 数 
为 count; 将 7 中 每 个 设施 有 服务 的 应 急 点 ex 都 放 入 对 应 的 
集合 Ri 中 ; 

d) 对 了 T 中 的 每 一 个 设施 f, 执行 如 下 操作 : 若 集 合 Ri 中 
的 任 一 应 急 点 ei 都 满足 min_q2(i, 7T) < bo， 则 把 设施 有 加 入 集 
合 M， 也 就 说 当 设 施 f 服 务 的 任 一 应 急 点 到 设施 集合 7 的 次 
小 距离 值 都 比 当 前 的 bo 小 ， 则 把 设施 加 入 集合 M; 若 集 合 
Ri 中 的 某 一 应 急 点 ei 满足 min_qz(i, 7) > bo， 则 把 设施 有 加 入 
集合 N， 也 就 说 当 设 施 f 服 务 的 某 一 应 急 点 到 设施 集合 7 的 
次 小 距离 值 比 当前 的 bo 大， 则 把 设施 f 加 入 集合 Ni 

e) 若 MM 不 为 空 ， 对 M 中 的 每 一 个 设施 fiw， 执行 如 下 操 
作 : 将 设施 fim 对 应 集合 中 的 Rm 中 的 所 有 应 急 点 eim 连接 到 
其 对 应 的 min_n2(im, 上 ,并 从 集合 T 中 删除 对 应 的 设施 fom， 
更 新 count， 更 新 每 个 应 急 点 到 集合 7 中 各 设施 的 最 短 与 次 
短 距离 ， 若 此 时 count<p， 则 此 时 求解 得 到 最 优 解 ，opt_flag 
= 1， 整 个 算法 结束 ; 若 count > Dp, 执 行 步 又 有 ; 若 M 为 空 ， 
执行 步 人 ; 

f) 更 新 初始 下 界 。 对 于 集合 中 每 一 个 集合 Ri 的 
max(mini_d2(i))， 将 这 些 最 大 值 max(mini_q2(i)) 中 的 最 小 值 令 
为 bs， 车 bi> bo， 则 更 新 下 界 为 bt; 若 bi< bo， 则 下 界 仍 为 bo 
= max(mini_d1(i))， 算 法 结束 。 
2.4 上 田子 算法 

事实 上 ， 任 何 一 个 可 行 解 对 应 的 目标 函数 值 均 能 作为 该 
问题 的 上 界 。 本 文 先 利用 前 面 所 介绍 的 数学 性 质 ， 将 问题 进 
行 降 阶 处 理 ， 之 后 执行 如 下 的 上 界 子 算法 求 出 上 界 : 

a) 对 于 所 有 ei:E€E EE， 按 照 其 mini_q2(i) 的 降序 排列 ， 排 序 
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c) 将 所 有 还 没有 被 服务 的 应 急 点 ex 与 集合 中 对 应 的 
设施 min_mi(is, 相连 接 ; 

d) 对 于 所 有 eiEE, 将 max(min_di(i, 加)=uo 作 为 上 界 ， 
令 max(min_di(is, UD)) 所 对 应 的 设施 设 为 fh; 

e) 若 集 合 信 中 存在 一 个 设施 ， 若 用 设施 来 蔡 换 开设 
的 设施 万 ， 即 执行 = UUf- 有 = VUf-fj， 之 后 如 果 对 
于 所 有 ei:E€ E，max(min_di(i, U)) < uo， 则 决定 用 设施 来 蔡 
换 开 设 的 设施 所 , 跳 到 步 虞 四 执行 ; 若 不 存在 这 样 的 设施 f， 
则 该 子 算法 结束 。 

3 ” 降 阶 回溯 算法 

本 算法 分 为 降 阶 算法 和 回溯 算法 两 个 部 分 ， 降 阶 算法 通 
过 结合 数学 性 质 进行 降 阶 ， 缩 小 问题 的 规模 ; 回溯 法 是 采用 
深度 优先 法 搜索 解 空间 。 回 渊 算法 对 设施 集合 FFs 中 的 每 一 
个 设施 分 为 以 下 2 种 情况 进行 处 理 : 9) 若 设施 万 开设 ，FF = 
FFIU 人 内 ，FF5s=FF5/ 内， 并 搜索 对 应 的 左 子 树 ; b) 若 设施 
矿 不 开设 ，Fro= FFoU {位 ，FFs = FFs / {f}， 并 搜索 对 应 的 
右 子 树 。 

3.1 降 阶 子 算法 

降 阶 子 算 法 步骤 如 下 : 

a) 初始 化 opt flag = 0, cur i=1, cur n=0, best gq= 
+oo; u_qg=0; 

b) 结合 前 文 的 数学 性 质 ， 若 某 个 设施 一 定 开设 ，F1 = 
UU {说 ，Fs= 了 ;1 {f}; 若 某 个 设施 一 定 不 开设 ，Fo= FoU 
{}，Fs = 二 Fs/ {f}。 若 执行 以 上 操作 后 集合 Fs 有 变化 ， 则 重 
新 执行 步 虞 a); 

c) 根据 前 文 的 上 界 子 算 法 对 问题 计算 上 界 , 此 时 开设 情 
况 赋 给 集合 Fpesr;， 对 应 的 目标 函数 值 赋 给 pest_9; 

d) 调用 下 一 节 介 绍 的 回溯 子 算法 。 
3.2 ”回溯 子 算法 

本 文 回溯 法 采用 深度 优先 法 搜索 问题 解 空间 ， 从 根 节点 
出 发 遍历 所 有 人 解 空间 ， 搜 索 至 任 一 节点 时 ， 判 断 该 节点 的 理 
论 上 界 pest_9 是 否 大 于 下 界 b_q， 若 满足 则 继续 疝 下 深度 优 
先 搜索 ， 否 则 进行 前 枝 ， 即 跳 过 该 节点 及 以 下 的 子 树 ， 逐 级 
向 上 回溯 。 
算法 执行 前 先 调用 前 面 介绍 的 降 阶 子 算 法 , 若 opt_1ag= 
1, 说 明 当 前 已 取得 最 优 解 , 算法 结束 ; 否则 执行 回溯 子 算法 。 
为 便于 描述 , 对 Fs 中 的 设施 重新 编号 , 初始 化 FFs5=Fs, FF 
= 人 ，FF0= 分 ， 再 调用 递归 回溯 子 程序 backtrack(1)。 

可 济 子 程序 backtrack(cur 让 描述 如 下 : 

a) 若 当 前 搜索 层 数 cur_i> | 上 | 或 者 IFFslF=0， 说 明 搜 索 到 
达 叶 子 节点 ， 此 时 若 | 而 |+|EEI<P， 则 得 到 一 个 可 行 解 。 若 
对 应 的 目标 函数 值 0 < best g， 则 更 新 最 优 解 ，best gq = 2， 
Fpest = Fi1U FF1; 结束 该 子 程 序 ; 

b) 检查 已 开设 设施 FiUFF1 是 否 已 覆盖 所 有 应 急 点 且 
EUFF < p， 若 满足 ， 则 达到 叶子 节点 ， 得 到 一 个 可 行 解 。 
若 对 应 的 目标 函数 值 O < best g， 则 更 新 最 优 解 ，opt_flag = 
1，best gq = OQ，Fbest= FiU FF1; 结束 该 子 程序 ; 

c) 在 假设 设施 jc; 开设 的 情况 下 ,计算 下 界 5b_g。 若 |FFI| 
<p-|mI 有 |FFs|>p-|fI-|FFI|， 同 时 bg < best gqg， 则 开设 
当前 设施 ， 执 行 FF1 = FFIU {fewr i，FFs=FFs/ {fcwri 调 


之 后 , 依 此 将 应 急 点 ei 连接 到 对 应 的 设施 mini_n1(i) 上 ， 逐 
开设 对 应 设施 ， 直 到 设 开设 个 数 为 p 个 或 者 所 有 应 急 点 都 被 
服务 为 止 ; 将 开设 的 设施 放 入 集合 U 中 ,将 还 没有 开设 的 的 
设施 放 入 集合 人 中 ; 

b) 若 所 有 的 应 急 点 全 部 被 集合 U 中 的 设施 服务 ， 则 此 
时 求解 得 到 最 优 解 ，opt_fiag = 1， 整 个 算法 结束 ;否则 继续 
执行 下 一 步 ; 


] backtrack(cur i+1) 进 入 左 子 树 搜 索 ; 若 上 述 条 件 不 成 立 ， 
则 跳 到 Step 6， 此 时 左 子 树 前 枝 ; 
d) 车 opt fiag=1, 说 明 当 前 算法 已 取得 最 优 解 ,， 结束 该 
子 程序 ， 否 则 继续 执行 该 子 程序 ; 

e) 算法 跳 到 搜索 树 上 一 层 前 , 恢复 FF1 和 Frs 到 步 又 c) 
的 初始 状态 ，FF = FF/ {feuri}， FFs= FFst+ {fouri}; 
f 在 假设 设施 fc ;不 开设 的 情况 下 ， 计 算 下 界 b5_q。 若 
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IFFi|<p-IFMAIFFs >p-|A|-|FFII， 同时 bgq<best qg， 则 
不 开设 当前 设施 ,执行 FF0= FFoU {fewi}, FFs= FFs/ {fe i}， 
调用 backtrack(cur i + 1) 进 入 右 子 树 搜 索 ， 若 上 述 条 件 不 成 
立 ， 则 结束 该 子 程序 ， 此 时 右 子 树 剪 枝 。 

g) 返回 搜索 树 上 一 层 前 , 恢复 FFo 和 FF; 到 步 肯 矣 的 初 
台 状 态 ，FFo = FFo/ {four i}， FFs= FFs+ {fowi}; 

算法 结束 后 ， 当 前 的 最 优 目标 函数 值 best gq 和 对 应 开设 
的 设施 Fyest 就 是 整个 问题 的 一 个 最 优 解 。 
3.3 ”算法 时 间 复 杂 度 分 析 

本 文 研究 的 问题 规模 即 设 施 个 数 为 n， 该 算法 在 搜索 空 
间 树 中 最 多 产生 2 个 叶子 节点 ,在 降 阶 子 算法 被 调用 后 ， 问 
题 规模 减少 为 |Fs|, 令 k=|Fs|， 因此 算法 在 最 坏 情况 下 的 时 间 
复杂 度 为 0(2”)。 

针对 应 急 选 址 问题 许多 学 者 提出 不 同 算法 ， 这 些 算 法 主 
要 分 为 精确 算法 、 近 似 算法 和 启发 式 算法 。 近 似 算法 与 启发 
式 算 法 的 优点 是 求解 速度 快 ， 不 足 之 处 在 于 一 般 无 法 得 到 问 
题 的 最 优 解 。 本 文 的 精确 算法 相对 于 一 般 的 精确 算法 而 言 ， 
通过 研究 其 数学 性 质 , 对 问题 进行 降 阶 ,缩小 问题 的 解 空间 ， 
调用 回溯 算法 前 对 搜索 树 进行 合理 剪 枝 ， 只 对 部 分 解 子 树 搜 
索 ， 加 快 了 求解 速度 。 


4 ”示例 分 析 及 实例 对 比 


4.1 示例 分 析 

为 了 更 清晰 的 说 明 本 文 算法 ， 下 面 给 出 一 个 示例 对 算法 
执行 的 整个 过 程 进行 说 明 。 

示例 1 如 图 1 所 示 ， 现 有 6 个 设施 j，5 个 应 急 点 鹿 
从 中 最 多 选 出 p=3 个 设施 , 使 得 最 大 需求 距离 最 小 , 表 1 给 
出 了 每 个 应 急 点 与 设施 之 间 的 距离 。 


图 1 应 急 点 与 设施 二 分 图 


Fig.1 Bipartite graph ofemergency points and facilities 
表 1 应 急 点 与 设施 之 间 的 距离 
Tab. 1 The distance between the emergency point and the facility 

J 

1 : 。 。 。 。 二 
1 了 ]3 ]4 js js 

11 3 5 7 6 8 10 

12 2 1 4 3 4 8 

js 6 5 6 7 8 9 

14 4 4 3 5 5 8 

ls 7 7 1 6 6 8 

i6 8 7 8 2 6 8 


对 示例 1 的 分 析 过 程 可 描述 如 下 : 
a) 根 据 性 质 4， 设 施 方 相 对 于 施 六 严格 劣势 ， 故 从 原 二 
分 图 中 删除 设施 js 并 且 删 除 边 {i1js，i2js，i3js，i4js， isjs，i6 


b) 根 据 性 质 4 或 者 性 质 5， 设 施 je 可 以 删除 ， 故 从 原 二 
分 图 中 删除 设施 je 并 且 删 除 边 {i1je，i2je，i3j6，i4j6, isj6, i6 


c) 根 据 性 质 6， 应 急 点 产 相 对 于 五 可 以 不 考虑 ， 故 将 应 
急 点 疡 删除 ， 从 原 二 分 图 中 删除 应 急 点 疡 并 且 删 除 边 { 广 ， 


i2j2, 12j3, i2j4, i2js, i2j6}; 
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中 根据 性 质 7， 应 急 点 i 相对 于 3 可 以 不 考虑 (因为 之 前 
设施 je 己 被 删除 )， 故 将 应 急 点 坟 4 删 除 ， 从 原 二 分 图 中 删除 应 
急 点 i， 并 且 删 除 边 {iaj1，iaj2，i4j3， i4j4, i4]J5s,， i4j6}; 

6) 经 过 降 阶 处 理 后 ， 更 新 问题 规模 如 图 2。 


三， 
> ,并 
CR NR oo Se 7 Oy 


3 6 


i i3 is i6 
图 2 问题 降 阶 处 理 后 的 二 分 图 

Fig.2 Bipartite graph of the problem after order reduction processing 
此 时 将 当前 剩余 应 急 点 与 剩余 设施 集合 中 距离 其 最 近 

的 设施 相连 ,计算 当前 开设 设施 数 得 4>p， 此 时 计算 出 下 界 
b=5。 
多) 执行 上 界 子 算法 ， 计 算出 上 界 x= 6。 


h) 执 行 回溯 算法 ， 其 执行 过 程 如 图 3 所 示 的 二 又 树 。 


图 3 解 空间 二 叉 树 
Fig.3 Solution space binary tree 

降 阶 后 的 解 空间 二 又 树 如 图 所 示 ， 通 过 上 述 算法 得 到 该 
问题 的 最 优 解 集 Prev= {j2，j3，j4}， 最 优 日 标 函 数值 best gq = 
5， 算 法 结束 。 
4.2 实验 结果 对 比 与 分 析 

本 文 实验 在 CPU 为 Core i7-7500 @2.70GHz， 操 作 系 统 
为 64 位 Windows10 的 环境 下 进行 。 采 用 Benchmark Library 
中 提供 部 分 基准 数据 进行 实验 ， 选 取 The P-median Problem 
实例 中 部 分 算 例 ，m=n=128， p=16， 该 数据 中 服务 距离 取 
值 为 (0, 4)， 采 用 Python 编程 实现 本 文 精确 算法 (BA) 和 蚁 群 
算法 (ACO)， 对 相同 数据 进行 求解 ， 并 对 两 种 算法 的 结果 进 
行 对 比分 析 , 实验 结果 记录 了 问题 求解 结果 和 程序 运行 时 间 ， 
如 表 2 所 示 。 


表 2 实验 结果 对 比 


Tab.2 Comparison of experimental results 


result Running time/s 
2 BA ACO BA ACO 
Pmp334 3 4 25.879467 1.212099 
Pmp434 4 4 25.254210 1.123443 
Pmp534 3 4 27.854497 1.032928 
Pmp634 3 4 25.309252 1.150121 
Pmp734 2 4 26.830491 1.230132 


从 示例 1 可 以 清楚 地 看 到 ， 利 用 数学 性 质 能 降低 问题 的 
规模 ， 加 快 算法 的 执行 速度 ， 从 上 述 实验 结果 可 得 ， 实 验 结 
果 与 理论 分 析 的 结论 一 致 : 蚁 群 算法 是 启发 式 算法 ， 计 算 速 
度 比 精确 算法 快 , 而 本 文 设计 的 算法 是 精确 算法 ,每 次 都 能 求 
得 最 优 解 ， 因 此 在 上 述 实验 中 本 文 算法 求 出 的 结果 稳定 性 要 
优 于 启发 式 算 法 类 型 的 蚁 群 算法 。 

针对 于 选 址 问题 设计 算法 的 国内 外 研究 一 般 分 为 启发 式 


\! 
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算法 、 近 似 算法 和 精确 算法 ， 本 文 提 出 的 降 阶 回溯 算法 是 精 
确 算法 ， 该 算法 利用 了 数学 性 质 、 上 界 及 下 界 来 加 快 算法 的 
运行 速度 ， 对 比 于 启发 式 算法 和 近似 算法 ， 本 文 设计 的 精确 
算法 保证 每 次 都 能 得 到 正确 稳定 的 全 局 最 优 解 ， 而 启发 式 算 
法 容易 陷入 局 部 最 优 ， 近 似 算法 也 只 能 以 一 定 的 相似 比 接近 
最 优 解 ， 一 般 无 法 得 到 问题 的 最 优 解 ， 该 问题 的 精确 算法 相 
关 研 究 内 容 较 少 , 与 一 般 的 精确 算法 相 比 最 终 的 解 是 一 样 的 ， 
但 是 一 般 精确 算法 都 是 对 所 有 节点 进行 分 析 ， 时 间 复 杂 度 比 
较 高 ， 本 文 利用 了 问题 的 数学 性 质 ， 结 合 设 计 的 上 下 界 子 算 
法 ， 能 够 批量 确定 某 些 开设 或 不 开设 的 设施 ， 降 低 问题 的 规 
模 ， 只 需 分 析 部 分 节点 ， 在 一 定 程度 上 降低 了 精确 算法 处 理 


= 


大 规模 问题 的 困难 ， 因 此 本 文 设计 的 算法 比 一 般 的 精确 算法 
速度 快 。 
5 ”结束 语 


本 文 研究 应 急 设施 选 址 问题 中 的 P- 中 心 模型 ， 针 对 启发 
式 算法 无 法 得 到 问题 最 优 解 和 精确 算法 处 理 大 规模 问题 具有 
较 高 时 间 复 杂 度 的 劣势 ， 本 文 设 计 了 一 种 能 快速 降低 问题 规 
模 同 时 能 得 到 精确 解 的 降 阶 回溯 精确 算法 ， 该 算法 的 设计 目 
的 在 于 保证 获得 全 局 最 优 解 的 情况 下 , 尽 可 能 减少 求解 时 间 。 
文章 通过 研究 P- 中 心 选 址 问题 的 数学 性 质 并 给 出 相应 证 明 ， 
这 些 数学 性 质 结合 所 设计 的 上 下 界 子 算法 ， 在 搜索 过 程 中 能 
批量 确定 某 些 开设 或 不 开设 的 设施 ， 在 回溯 过 程 中 进行 大 量 
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剪 枝 ， 从 而 降低 问题 求解 规模 ， 提 高 了 算法 的 执行 速度 ， 同 
时 ， 这 些 数 学 性 质 也 可 以 应 用 于 其 他 算法 中 。 通 过 理论 分 析 
及 示例 1 可 以 看 出 , 采用 本 文 的 算法 求解 P- 中 心 选 址 问题 时 
不 仅 能 求 出 精确 解 ， 同 时 利用 数学 性 质 和 上 下 子 算法 能 降低 
问题 的 规模 ， 加 快 问题 求解 速度 。 
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